


Steven J. Wagstaff,’ Michael J. Bayly,’ 
Philip J. Garnock-Jones,' 
and Dirk C. Albach?^ 


CLASSIFICATION, ORIGIN, 
AND DIVERSIFICATION OF 
THE NEW ZEALAND HEBES 
(SCROPHULARIACEAE)! 





ABSTRACT 


The New Zealand hebes (Scrophulariaceae) are members of a large Southern Hemisphere clade nested within Ve- 
ronica. Analysis of ITS and rbcL sequences suggests that the New Zealand species are derived from a single common 
ancestor that arrived via long-distance dispersal. After the establishment of this initial founder population in New 
Zealand, the hebes have undergone at least two major episodes of diversification, giving rise to six clades. The great 
degree of morphological diversity in the New Zealand hebes contrasts with a corresponding low level of sequence 
divergence. New Zealand was a source of new emigrants to other regions in the South Pacific that were preadapted to 
high mountains or forest margins. Our results suggest that two instances of long-distance dispersal from New Zealand 
to South America, at least one instance from New Zealand to Australia, and one instance from New Zealand to New 
Guinea have occurred relatively recently. Shorter hops to the Chatham Islands and the subantarctic islands are also 


supported by the sequence data. 
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Long-distance dispersal has a profound influence 
on the evolution of insular floras (Carlquist, 1974), 
and there is substantial evidence suggesting that it 
occurs. relatively frequently (Godley, 1967; Pole, 
1994). One of the most remarkable examples of dis- 
persal followed by adaptive evolution on islands is 
the New Zealand hebes (Scrophulariaceae). Wags- 
taff and Garnock-Jones (1998, 2000) suggested that 
the New Zealand hebes are the descendants of a 
small founder population that may have been de- 
rived from a single seed. They proposed that com- 
bined influences of inbreeding, genetic drift, and 
strong selection acting upon small populations have 
probably played a major role in the rapid diversi- 
fication of the group. 

The hebes are one of the largest and most eco- 
logically diverse plant groups in New Zealand, in- 
cluding over 120 species, with outlier populations 
in eastern Australia, Tasmania, New Guinea, Rapa 
Island, and South America. They range from al- 
pine cushion-forming plants (Fig. 5D) to lowland 
woody shrubs or small trees (Fig. ЭМ, К, 5), and 


in New Zealand are conspicuous elements in most 
terrestrial ecosystems except forests and wetlands. 
Species such as Hebe armstrongit, H. cupressoides, 
and H. speciosa have patchy or localized distri- 
butions and are considered rare or endangered; 
about 70% of the species are confined to small 
regions within New Zealand. 

The New Zealand hebes were formerly included 
in a broadly defined circumscription of the genus 
Veronica (Wettstein, 1891; Cheeseman, 1925), but 
recent flora and taxonomic treatments (Ashwin & 
Moore in Allan, 1961; Garnock-Jones, 1993а, b; 
Heads, 1994a, b) recognize less inclusive groups 
(see Table 1), usually accepting four genera in New 
Zealand: Chionohebe, Hebe, Heliohebe, and Para- 
hebe (Garnock-Jones, 1993a, b). Heads (1987) de- 
scribed an additional genus, Leonohebe. Although 
we do not accept his wide circumscription of that 
genus, the name Leonohebe could be applied to a 
small clade of four or five species that is supported 
by the analyses of Wagstaff and Garnock-Jones 
(1998, 2000). 
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Classifications of New Zealand hebes. The New Zealand species were placed in three sections of Wettstein's 


(1891) treatment of Veronica, and Cheeseman (1925) placed them in two divisions of Veronica. 


Ashwin & Moore 


Wettstein (1891) Cheeseman (1925) 


Veronica Veronica 
sect. Pygmea Division Pygmea Pygmea 
sect. Hebe 

Division Hebe Hebe 


"Semiflagriformes" 


"Connatae" 


"Flagriformes" 


"Buxifoliatae" 


“Subdistichae” 
"Subcarnosae" 


in Allan (1961) 


Heads (1987, 1994b) Garnock-Jones (1993a, b) 


Chionohebe Chionohebe Including 
Parahebe “Group B”) 
Leonohebe 


sect. Densifoliae 

Hebe 
sect. Leonohebe "Semiflagriformes" 
sect. Connatae "Connatae" 
sect. Apiti 
sect. Flagriformes — "Flagriformes" 
sect. Aromaticae 
sect. Salicornioides 
sect. Buxifoliatae — "Buxifoliatae" 
Hebe 
sect. Subdistichae sect. Subdistichae 
sect. Glaucae sect. Glaucae 


sect. Hebe sect. Hebe 


sect. Chamaedrys 


sect. Labiatoides 

sect. Paederota 

sect. Paederotoides 
sect. Pseudolysimachia 
sect. Veronicastrum 
sect. Omphalospora 
sect. Beccabunga 


"Apertae" 
"Occlusae" 


"Grandiflorae" 
"Paniculatae" 
Parahebe 
"Group A, B, C" 


Division Euveronica 


ser. Hebe 
ser. Occlusae 
Parahebe 


sect. Paniculatae 


ser. Hebe 
ser. Occlusae 


*Grandiflorae" 
Heliohebe 
Parahebe 
"Group А, С“ 
Derwentia 





! In the Flora of New Zealand Volume | (Allan, 1961), M. B. Ashwin prepared the treatment of Parahebe, Рудтеа, 
and the informal grouping "Flagriformes" of Hebe. The remainder of the Hebe treatment, including the informal synopsis, 


was prepared by L. B. Moore. 


This research contributes to ongoing efforts to 
create a phylogenetic classification of Scrophulari- 
aceae. Olmstead and Reeves (1995) and Olmstead 
et al. (2001) showed that the Scrophulariaceae, as 
traditionally cireumscribed, are not monophyletic. 
They identify clades from a dismembered Зсгорћи- 
lariaceae s.l. that could merit formal recognition. In 
their studies Veronica was nested within a large 
clade they called the Antirrhinaceae nom. cons. 
prop. (Reveal et al., 1999). This large clade was 
recognized by Olmstead and Reeves (1995) and in- 
cludes part or all of Bentham’s (1876) tribes Digi- 
taleae, Antirrhineae, Cheloneae, and Gratioleae. 
the small tribe Angelonieae, and the small families 
Callitrichaceae, Globulariaceae (excluding Selagi- 
naceae), Hippuridaceae, and Plantaginaceae. 

The aim of this research is to identify well-sup- 


ported monophyletic groups among the New Zea- 
land hebes, to improve their classification, infer 
their origin, and explore underlying processes of 
diversification. We propose that diversification in 
the group reflects transoceanic dispersal and adap- 
tive radiation. The hebes have successfully exploit- 
ed a diversity of ecological niches that were prob- 
ably created during the recent uplift and glaciation 
of the mountains of New Zealand. 


MATERIALS AND METHODS 


Our sampling strategy capitalized on the unique 
characteristics of rbcL and ITS sequences. The 
plastid encoded gene rbcL has relatively few vari- 
able sites, which allowed sequence comparisons 
among distantly related outgroups, and placement 
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of the hebes within Scrophulariaceae. It is also use- 
ful in that a large number of published rbcL se- 
quences are available for comparison (see Chase et 
al., 1993; Küllersjó et al.. 1998; Olmstead et al., 
2001, and references therein). Finally. Albert et al. 
(1994) and Bremer and Gustafsson (1997) suggest- 
ed that the gene rbcL approaches clock-like behav- 
ior in its evolution, and hence the amount of se- 
quence divergence could be used to estimate 
divergence times. Ву comparison, the nuclear en- 
coded ITS-region has many more variable sites than 
rbcL, which provides more informative characters 
to resolve relationships at lower taxonomic levels 
(Baldwin et al., 1995). 


STUDY GROUP 


The rbcL study group consisted of 33 species 
including 12 of the New Zealand hebes with at least 
| representative from each of the currently recog- 
nized genera, 5 species of Veronica, and 1 species 
of Veronicastrum. Nineteen rbcL sequences were 
newly published herein along with 13 published 
sequences of Antirrhinaceae from Olmstead et al. 
(2001), and Nicotiana tabacum (Solanaceae) was 
designated as the outgroup (Lin et al., 1986). Seven 
sequences were considered redundant; even though 
they were not identical, the resolution of missing 
data could potentially make them identical. We 
therefore excluded Derwentia derwentiana and D. 
perfoliata, and the New Zealand accessions of Hebe 
elliptica and H. salicifolia, from subsequent analy- 
ses. Nineteen of the 37 species included in the rbcL 
analysis were also included in the ITS survey. 

The ITS study group included 78 sequences, 19 
of which were newly published. Among these are 
58 representatives of Chionohebe, Derwentia, He- 
liohebe, and Parahebe, including conspecific acces- 
sions of Chionohebe ciliolata and C. densifolia from 
Australia and New Zealand, Parahebe lithophila 
from Australia, and P. randewateri from New Guin- 
ea. Thirty-five species of Hebe were also included, 
and among these were: at least one representative 
from each of Moore's (in Allan, 1961) informal 
groups; H. formosa from Tasmania: H. benthamii 
from the New Zealand subantarctic islands; H. bar- 
keri, H. chathamica, and H. dieffenbachii from the 
Chatham Islands (east of the main islands of New 
Zealand); accessions of H. elliptica from both New 
Zealand and the Falkland Islands; and H. salici- 
folia from both New Zealand and Chile. Pseudoly- 
simachion, Veronica, Veronicastrum, and Wulfenia 
emerged as potential sister groups of the hebes in 
the analysis of Hong (1984) and Albach and Chase 
(2001); therefore a total of 19 species representing 


these. genera. were included in our analysis. The 
Asiatic species Veronicastrum sibiricum was desig- 
nated as the outgroup for the analysis of ITS se- 
quences. 

Voucher specimens are listed in Appendix 1, 
along with collection information, literature cita- 
tions, and GenBank (<http://www.nebi.nim.nih. 
gov) accession numbers. The complete data sets 
are available upon request from the first author, and 
they were deposited in TreeBASE (http://www. 
herbaria.harvard.edu/treebase-). The study acces- 
sion number is 5623. and the matrix accession 
numbers are M961 (rbcL) and M962 (ITS). 


DNA EXTRACTION, AMPLIFICATION, AND 
SEQUENCING 


Total DNA was extracted from either fresh leaves 
or leaf fragments dried with silica gel using a mod- 
ification of the hot CTAB method of Doyle and 
Doyle (1987). The cpDNA gene rbcL and the 
nrDNA ITS-region [the 3' end of the 185 rDNA 
gene; internal transcribed spacer —1 (ITS-1): the 
5.88 rDNA gene: internal transcribed spacer ~2 
(ITS-2); and the 5' end of the 285 rDNA gene] were 
amplified by PCR. Primer sequences and our am- 
plification and sequencing techniques follow Olm- 
stead et al. (1992) Гог rbcL, and Wagstaff and Gar- 
nock-Jones (1998) for the ITS-region. Excess 
primers and unincorporated nucleotides were re- 
moved from the PCR products by spin column cen- 
trifugation (QlAquick PCR purification kit, QIA- 
GEN Inc.). The purified DNA samples were then 
labeled with Big Dye terminators (PE Applied Bio- 
systems, The Perkin-Elmer Corp.). Both the forward 
and reverse DNA strands were sequenced by the 
Waikato University DNA Sequencing Facility. Con- 
tig editing and assembly was accomplished using 
Sequencher version 3.0 (Gene Codes Corp.). 


SEQUENCE ALIGNMENT 


The sequence alignment for the ITS-region was 
facilitated by ClustalX (Thompson et al., 1997). A 
gap penalty setting of 75 and a gap extension pen- 
айу of 6.6 were initially used to identify and. po- 
sition large gaps in the sequence data; then low- 
scoring segments were realigned using a gap 
penalty setting of 15 and a gap extension penalty 
of 6.6 with the removing new gaps option turned 
on. These settings opened and positioned small 
gaps. The final alignment was inspected and minor 


revisions were made manually. 
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DATA ANALYSIS 


The phylogenetic analyses were accomplished 
using PAUP* version 4.0d65 (Swofford, 1998). The 
analyses were conducted using the PAUP* settings 
random addition sequence with 100 replicates. 
TBR branch swapping. mulpars in effect. and 
steepest descent. The characters were all unordered 
and weighted equally, and gaps were treated as 
missing data. 

Support for the inferred clades is given by jack- 
knife percentages (Farris et al., 1996) and by the 
number of synapomorphies for each group. Jack- 
knife analysis investigates the structure, or phylo- 
genetic signal, in a matrix without permutation, but 
excludes an assigned fraction of characters, here 
set to 30%. The jackknife searches were performed 
with 1000 replications excluding uninformative 
sites, maxtrees = 10 for each replication; the start- 
ing trees were obtained by random addition with 
one replication for each jackknife replication, TBR 
branch-swapping, and mulpars in effect. 

The relationship between sequence divergence 
and time for the gene rbcL was discussed by Albert 
et al. (1994) and Bremer and Gustafsson (1997) and 


was calculated using the equation: 


substitution rate = patristic distance(Dp)/number 


of nucleotides/inferred time since cladogenesis. 


RESULTS 


The rbcL sequences were 1402 nucleotides in 
length (positions 27-1428 in tobacco). Among the 
1402 sites included in the rbcL matrix, 1103 sites 
were invariant; 144 were parsimony-uninformative, 
and 155 characters were parsimony-informative. 
Missing data accounted for 4.7% of the matrix. 
Most of the missing data fell in a region at the 5' 
end of rbcL, upstream from the conserved EcoRV 
restriction. site used in cloning some sequences 
(Olmstead et al.. 1992), and in a region downstream 
from position 1325, the location of a PCR primer 
site used to amplify some sequences. Most changes 
(calculated across the maximum parsimony tree 
shown in Fig. 2) occurred in the third codon posi- 
tion (309): substantially fewer changes occurred in 
the first (116) and second (17) codon. positions. 
Changes in the third codon position are generally 
synonymous, and hence are more likely to evolve 
in a clock-like manner. 

The analysis of rbcL sequences recovered 9408 
maximum parsimony trees distributed in a single 
island of 529 steps (consistency index = 0.52 ex- 
cluding uninformative characters, retention. index 
= 0.73): a strict consensus tree is shown in Figure 


l. and one of the maximum parsimony trees 15 
shown in Figure 2. Most members of the Antirrhin- 
aceae form a basal grade in our analysis with Chio- 
nohebe, Derwentia, Hebe, Heliohebe, Parahebe, and 
Veronica forming a clade that receives 100% jack- 
knife support. V. anagallis-aquatica is sister to a 
largely Australasian clade (98% jackknife support) 
that includes the New Zealand hebes along with 
Derwentia, Veronica arguta, and V. persica. Rela- 
tionships within this clade are poorly resolved 
(Figs. 1, 2). 

The mean absolute distance and standard devi- 
ation from Veronicastrum sibiricum to the Chionoh- 
ebe, Derwentia, Hebe, Heliohebe, Parahebe, and Ve- 
is 29.1 + 7.7 (see Fig. 2). 


Veronicastrum is reported in the fossil record (Tiff- 


ronica terminals 
пеу, 1985) from the mid Miocene some 15 million 
vears before present (mybp). The substitution rate 
in the Veronicastrum lineage was estimated by di- 
viding 29.1 + 7.7/15 = 1.9 + 0.5 substitutions for 
the entire gene rbcL. The mean distance from the 
terminals to the ancestral node of the Australasian 
species (including V. persica, which is Eurasian) is 
18.8 changes, which corresponds to an upper Mio- 
cene divergence estimate of about 9.9 mybp. The 
mean distance from the terminals to the ancestral 
node of the H. salicifolia lineage is 7.4 changes. 
which corresponds to a Pliocene divergence esti- 
mate of about 3.9 mybp (Fig. 2). 

The aligned ITS matrix was 695 nucleotides in 
length with gaps created to account for insertions 
and deletions, among which 364 sites were con- 
stant, 109 were parsimony-uninformative, and 222 
were potentially parsimony-informative. Missing 
data accounted for 1.296 and gaps accounted. for 
9.6% of the ITS data matrix. The 5.85 gene was 
uniformly 165 nucleotides; ITS-1 varied between 
175 and 224 nucleotides and ITS-2 between 203 
and 215 nucleotides. Most of the variation in the 
ITS region was observed in ITS-1 and ITS-2. The 
5.88 gene was more conserved. Conserved motifs 
identified by Liu and Shardl (1994) and Hershkov- 
itz and Zimmer (1996) were identified in the ITS- 
| and ITS-2 sequences in our survey. 

Thirty-eight insertions and deletions (indels) 
were inferred in the ITS-1 and ITS-2 spacer regions 
(Table 2). Mostly these were relatively small. rang- 
ing from 1 to З bp, but two large insertions of 27 
and 44 bp and three deletions of 7 or 8 bp were 
also inferred. Most of the indels were unique to a 
single sample, but 13 were shared by two or more 
species, sometimes uniting groups supported by 
substitutions in the sequence data alone. Species 
of Derwentia are characterized by a one-base de- 
letion. Species of Heliohebe are characterized by a 
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Hebe elliptica 

Hebe benthamii 

Hebe odora 

Hebe salicifolia 

Hebe salicornioides 
Parahebe vandewateri 
Hebe macrantha 

Hebe cupressoides 
Heliohebe raoulii 
Chionohebe densifolia 


Hebe cheesemanii 


New Zealand hebes 


Parahebe catarractae | 
Derwentia nivea 
Veronica arguta 
Hebe formosa 
Veronica catenata 
Veronica persica 
Veronica anagallis-aquatica 
Veronica officinalis 
Digitalis purpurea 
Plantago lanceolata 
Veronicastrum sibiricum 
Callitriche heterophylla 
Hippuris vulgaris 
Chelone obliqua 
Collinsia grandiflora 
Antirrhinum majus 
Globularia cordifolia 
Gratiola pilosa 
Amphianthus pusillus 
Bacopa caroliniana 
Angelonia pubescens 


Nicotiana tabacum 


Figure 1. Strict consensus of 9408 minimal length trees produced by parsimony analysis of rbcL sequences. This 
tree shows the placement of New Zealand hebes within the Antirrhinaceae sensu Olmstead et al. (2001). using Nicotiana 
tabacum as an outgroup. Jackknife values > 50% are given above each node. 
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Hebe benthamii 
Hebe salicornioides 
Hebe odora 
“ Hebe elliptica 
3.9 Hebe salicifolia 
mybp Derwentia nivea 
Veronica arguta 
Parahebe catarractae 
Parahebe vandewateri 
Hebe cupressoides 
Hebe macrantha 
Chionohebe densifolia 
Veronica catenata 
Veronica persica 
Hebe formosa 
v Hebe cheesemanii 
Heliohebe raoulii 


9.9 
mybp 


Veronica anagallis-aquatica 
Veronica officinalis 
15 Veronicastrum sibiricum 
mybp Digitalis purpurea 
Plantago lanceolata 
Callitriche heterophylla 
Hippuris vulgaris 
Globularia cordifolia 


Chelone obliqua 


Collinsia grandiflora 
Antirrhinum majus 


Gratiola pilosa 
Amphianthus pusillus 
Bacopa caroliniana 


Angelonia pubescens 


Nicotiana tabacum 
— 5 changes 


Figure 2. One of the maximum parsimony trees recovered from a parsimony analysis of rbcL sequences. The branch 
lengths are proportional to the number of changes along each branch. See scale at bottom. Fossils of Veronicastrum 
sibiricum are reported from the mid Miocene about 15 mybp. Divergence estimates are provided at the ancestral node 
of Australasian species (V. persica is Eurasian) and the ancestral node of the Hebe clade. 
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Table 2. Insertions and deletions inferred from ITS sequence comparison. 


— — —— — — — ——— —— —— — — — —— — 
— ш — — — — TT — — — — — — — —— = 


Insertion/ 
Taxon deletion Size Position Sequence 
Veronicastrum sibiricum Deletion 1 Ч Uc 
Veronicastrum sibiricum, Wulfenia Insertion | 54 e 
carinthiaca 
Veronica glandulosa Deletion | 61 Ug 
Hebe vernicosa Insertion 3 68—79 gla 
Veronica calycina Insertion | 17 ( 
Veronica chamaedrys Deletion | 78 c/g 
Veronicastrum sibiricum Deletion | 89 с 
Wulfenia carinthiaca Insertion 44 90-134 aatctaggtgtgcaageccetttgttgagag- 
lecegegectgete 
Pseudolysimachion spicata, Veronica Insertion 27 108-134 gactaglegagtgegecgetetegete 
anagallis-aquatica, V. bellidioides, V. 
glandulosa, V. glauca, V. officinalis, V. 
serpyllifolia, V. urticifolia 
Veronica bellidioides, V. officinalis, V. Insertion l 162 a 
urticifolia 
Derwentia derwentiana, D. nivea, D. Insertion l 207 a 
perfoliata, Parahebe lithophila, 
Veronica arguta 
Veronica macrostachya Deletion 3 208-209 ес 
Hebe salicifolia Insertion | 217 с 
Veronica bellidioides, V. glandulosa, Insertion | 234 cg 
V. officinalis, V. urticifolia, Wulfenia 
carinthiaca 
Veronica persica Deletion | 248 а 
Veronica chamaedrys, V. anagallis- Insertion 1 200 а 
aquatica 
Pseudolysimachion spicata, Veronica Insertion | 442 a 
chamaedrys, V. macrostachya, V. 
oltensis 
Veronica urticifolia Insertion | 154 e 
Veronica anagallis-aquatica Insertion | 107 l 
Pseudolysimachion spicata, Veronica Deletion l 474 g 
anagallis-aquatica, V. bellidioides, V. 
fruticulosa, V. glandulosa, V. glauca, 
V. officinalis, V. saturejoides, V. 
serpyllifolia, V. urticifolia, 
Veronicastrum sibiricum, Wulfenia 
carinthiaca 
Parahebe canescens Insertion | 475 g 
Veronica bellidioides Insertion 2 501 се 
Wulfenia carinthiaca Deletion l 503 t 
Heliohebe hulkeana, H. lavaudiana, H. Insertion l 508 І 
raoulii 
Veronica persica Deletion | 513 [ 
Hebe elliptica, H. elliptica. Parahebe Deletion 8 583-590 tclegtge 
birleyi, P. brevistylis, P. decora, P. lyallii, 
Р. spathulata, P. vandewateri 
Parahebe planopetiolata Deletion 7 590—596 саїсїсс 
Parahebe canescens Deletion 8 591-598 alctecge 
Veronica persica Deletion 3 592-594. (са 
Heliohebe hulkeana, Н. lavaudiana, Н. Deletion | 606 g 
raoulii 
Parahebe vandewateri Deletion l 609 a 
Veronica chamaedrys Deletion 2 609—610 ag 


Eee 
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Table 2. Continued. 
Insertion/ 
Taxon deletion Size Position Sequence 

Parahebe canescens Insertion 3 615-616 cat 

Veronica chamaedrys Deletion 2 620-621 te 

Pseudolysimachion spicata Deletion 2 625—626 ac 

Veronica chamaedrys, V. serpyllifolia Insertion 1 631 a 

Veronica austriaca, V. oltensis Insertion 2 638—639 gc 

Veronica arguta Insertion l 640 c 





one-base insertion and a one-base deletion. Both 
accessions of Hebe elliptica and nine species of 
Parahebe have an eight-base deletion that appears 
to have evolved independently at least three times 
(Table 2). This deletion is lacking in Parahebe lin- 
ifolia and in P. catarractae subsp. catarractae and 
subspecies martinii. The South American accession 
of Hebe salicifolia has a unique one-base insertion 
that is lacking in the accession of Hebe salicifolia 
from New Zealand. 

The ITS sequences in our study were evolving at 
a faster rate than rbcL. The average rate of change 
per variable site for rbcL was 1.7 (tree length of 
532/number of variable sites 144 + 155). The 
aligned ITS sequences were shorter than rbcL, 
there were more variable sites (331), and the av- 
erage rate of change per variable site was 3.6 for 
the ITS region. 

Parsimony analysis of the ITS-region recovered 
6931 maximum parsimony trees distributed in at 
least two islands of 1213 steps (consistency index 
— 0.41 excluding uninformative characters, reten- 
tion index — 0.73); a strict consensus tree is shown 
in Figure 3 and one of the maximum parsimony 
trees in Figure 4. Based upon the results from anal- 
ysis of rbcL sequences (Figs. 1, 2). Veronicastrum 
sibiricum was designated as the outgroup. The ear- 
liest divergence within the ingroup is between Wul- 
fenia carinthiaca and all other taxa. The Northern 
Hemisphere species of Veronica are found in five 
clades that form a grade basal to a Southern Hemi- 
sphere clade comprising the New Zealand hebes 
and their relatives. A heterogeneous Australian 
clade comprised of Derwentia, Hebe formosa. Par- 
ahebe lithophila, Veronica arguta, and V. calycina 
(8296 jackknife; 18 synapomorphies) is sister to the 
New Zealand hebes, though there is relatively little 
support for this relationship (70% jackknife: 10 
synapomorphies) (Figs. 3, 4). 

Six well-supported clades are identified among 
the New Zealand hebes, but the relationships 
among these clades are unclear (Figs. 3. 4). The 


first is a clade that comprises Leonohebe s. str. in- 
cluding a well-supported group, Hebe tetrasticha, Н. 
cheesemanii, and Н. ciliolata (99% jackknife: 6 
synapomorphies) with Hebe cupressoides weakly 
supported as their sister (61% jackknife: 3 syna- 
pomorphies). The Chionohebe A clade (98% jack- 
knife; 5 synapomorphies) includes Parahebe plan- 
opetiolata and the cushion-forming species of 
Chionohebe with both the New Zealand and Aus- 
tralian accessions of C. ciliolata. The Chionohebe 
B clade consists of Parahebe trifida and both the 
New Zealand and Australian accessions of Chio- 
nohebe densifolia (91% jackknife: 4 synapomor- 
phies). The fourth clade includes 6 species of Par- 
ahebe (9396 jackknife; 7 synapomorphies) and 
accommodates the informal “Groups A & С“ of 
Ashwin (in Allan, 1961; Table 1) and P. spathulata. 
The fifth clade includes all the species of Heliohebe 
in our analysis (10096 jackknife: 13 synapomor- 
phies). The sixth includes the remaining species of 
Hebe with both New Zealand and South American 
accessions of Hebe salicifolia and Hebe elliptica 
(100% jackknife; 11 synapomorphies). Hebe ma- 
crantha is weakly supported as the sister to the rest 
of this clade (50% jackknife; 4 synapomorphies) 
(Figs. 3, 4). 


DISCUSSION 


Large. unwieldy genera with a cosmopolitan dis- 
tribution such as Veronica pose among the most dif- 
ficult taxonomic problems for plant systematists. 
whose opinions are often strongly held. One of the 
most vexing of these problems is the inconsistent 
means by which taxonomists define generic bound- 
aries and the recognition of rank within a hierar- 
chical classification scheme. Recent taxonomic 
treatments in the Southern Hemisphere have fa- 
vored narrow circumscriptions, and several new 
genera have been segregated from Veronica (see Ta- 
ble 1), whereas taxonomists in Europe and North 


America have traditionally embraced a broad ge- 
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Figure 3. Strict consensus of 6931 minimal length trees produced by parsimony analysis of the entire ITS-region. 
Notable clades are identified with brackets. Jackknife values > 50% are given above each node. 'Parahebe catar- 
ractae subsp. martinii, *Parahebe catarractae subsp. catarractae, ‘Hebe hectorii subsp. subsimilis, ‘Hebe hectorii 


subsp. hectorti. 


neric definition of Veronica, e.g., Wettstein (1891). 
This discrepancy of opinion contributes to taxo- 
nomic ambiguity and instability. 

Our results support those of Albach and Chase 
(2001), implying that the genus Veronica is at best 
paraphyletic by exclusion of the Southern Hemi- 
sphere genera Chionohebe, Derwentia, Hebe, He- 
liohebe, Leonohebe, and Parahebe (Figs. 1, 3). as 


well as the Eurasian genera Paederota and Pseu- 
dolysimachion, and the North American genera 
Synthyris and Besseya. One possible solution is to 
lump them all in a broad circumscription of Ve- 
ronica. This move. however, would create a cas- 
cade of nomenclatural changes requiring the rec- 
ognition of many new combinations and the 
adoption of old combinations within Veronica. We 
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Figure 4. Опе of the maximum parsimony trees recovered from a parsimony analysis of the ITS-region. The New 
Zealand hebes comprise a well-supported monophyletic group with outliers on the offshore islands and in Australia, 
New Guinea, and South America. Dispersal away from the main islands of New Zealand is inferred in nine species. 
Branch lengths are proportional to the number of changes along each branch. See scale at bottom. 'Parahebe catarractae 
subsp. martinii, "Parahebe catarractae subsp. catarractae, *Hebe hectorii subsp. subsimilis, *Hebe hectorii subsp. hectorii. 


accept that retaining a paraphyletic Veronica ob- 
scures phylogenetic relationships: however, an al- 
ternative approach is to recognize smaller, less in- 
clusive clades as generic segregates of Veronica. 
This approach was adopted by Hong (1984). Here 
we identify major clades supported by the DNA 


sequences and describe patterns of diversification 
in the New Zealand hebes. 


MAJOR CLADES OF NEW ZEALAND HEBES 


A heterogeneous clade composed of Derwentia, 
Hebe formosa, Parahebe lithophila, Veronica arguta. 
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Figure 5. Plate illustrating some of the morphological diversity in the New Zealand hebes. А. Flowering shoot 


of Derwentia perfoliata with toothed leaves obscure in this picture. — В. Flowers of Parahebe catarractae. —C. Panic- 


ulate and terminal inflorescence of Heliohebe raoulii subsp. maccaskillii, P. J. Garnock-Jones 2123. —D. Cushion habit 
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and V. calycina is sister to the New Zealand hebes 
(Figs. 3, 4). The genus Derwentia includes nine 
currently accepted species that are endemic to 
Australia, from southeastern Queensland to Tas- 
mania and west to Kangaroo Island in South Aus- 
tralia, where they are found mostly in tableland or 
cool temperate regions (Briggs & Ehrendorfer, 
1992). Toothed leaves [only on the lower branches 
of H. formosa and often obscure in D. perfoliata 
(Fig. 5A)]. and a dense ring of hairs in the corolla 
throat (glabrous in H. formosa) are possible synap- 
omorphies that unite the clade. Derwentia and Hebe 
formosa also have similar growth forms; their new 
shoots are initiated at the base of the plant, over- 
topping older, short-lived branches. Hebe formosa 
is distinguished from Derwentia by the occurrence 
of a one-base insertion (Table 2). The chromosome 
number of Hebe formosa is n = 21, this possibly 
being the ancestral state within the Derwentia 
clade; chromosome numbers of n = 19 or 20 are 
published for Derwentia (Briggs & Ehrendorfer, 
1992). 

We identify six major clades within the New Zea- 
land hebes that are supported by the sequence 
data, which we refer to as the Leonohebe clade, the 
Chionohebe A clade, the Chionohebe B clade, the 
Parahebe clade, the Heliohebe clade, and the Hebe 
clade (Fig. 3). ITS sequences provide strong sup- 
port for the Leonohebe clade (Figs. 3, 4), comprising 
Hebe cheesemanii (Fig. 5K), H. ciliolata (Fig. 5J). 
and Н. tetrasticha. Hebe cupressoides (Fig. 5L) is 
weakly supported as sister to this clade. These spe- 
cies are endemic to the South Island of New Zea- 
land. They have traditionally been included in 


Hebe, but here, as in previous analyses (Wagstaff & 
Garnock-Jones, 1998, 2000), are far removed from 
the other species that Moore (in Allan, 1961) in- 
cluded in that genus. Hebe cupressoides and mem- 
bers of the Leonohebe clade lack the dorsal capsule 
compression typically found in members of Hebe, 
and a potential synapomorphy for the group is pos- 
session of a distinctive cupressoid growth habit, 
which has apparently evolved independently from 
that in the whipcord hebes [e.g.. H. annulata, Н. 
armstrongii, H. hectorii, H. lycopodioides, H. pro- 
pinqua, and Н. salicornioides, which have a similar 
growth form (Figs. 5T-W)|. Hebe cheesemanii, H. 
ciliolata, and H. tetrasticha, along with H. tumida, 
comprise the informal group “Semiflagriformes” of 
Moore (in Allan, 1961). These similar species are 
subshrubs that occur in rocky areas at high altitude 
and are characterized by possession of lateral in- 
florescences (Fig. 5K), dioecious breeding system, 
and sour-scented flowers. We refer to this group as 
the Leonohebe clade, because Leonohebe ciliolata 
(H. ciliolata) was designated as the nomenclatural 
type for that genus by Heads (1987). If the group, 
including H. cupressoides, is treated as generically 
distinct from Hebe (and other genera of the Hebe 
complex; see Garnock-Jones, 1993a), then the use 
of the name Leonohebe seems warranted, though 
this is a much more restricted use of the name than 
the clearly polyphyletic circumscription originally 
employed by Heads (1987). 

The Chionohebe A clade encompasses the cushion- 
forming species of Chionohebe. All of these species 
occur in the South Island where they are high-alpine 
plants of rock and scree. The cushion-forming species 
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of Chionohebe thomsonii, male plant, Eyre Mis., South Island, P. J. Garnock-Jones 1906. —E. Lateral view of 
erect tubular flower of Chionohebe pulvinaris, from Takitimu Range, South Island. F, G. Hebe macrantha var. 
brachyphylla. —F. Habit of a plant ca. 15 ст tall from same population as M. J. Bayly 560, Mt. Arthur, South 
Island. —G. Apex of vegetative shoot from Mt. Arthur, WELT 82554. H, 1. Hebe epacridea. —H. Habit of a plant 
ca. 10 em tall, same population as M. J. Bayly 795, Mt. St. Patrick, South Island. —1. Branchlet of plant from unknown 
locality. —J. Hebe ciliolata, branchlet showing ciliolate leaf margins, plant from Mt. Arthur, WELT 82556. —K. Hebe 
cheesemanii, male plant, from same population as M. J. Bayly 756-757. Black Birch Ra., South Island. —L. Hebe 
cupressoides, branchlet, cultivated Atawhai, Nelson, South Island, WELT 82553. M, N. Hebe parviflora. —M. Habit of a 
plant са. 2.5 m tall from same population as P. J. Garnock-Jones 2258, Hauhangaroa Ra., North Island. —N. Inflorescence, 
P. J. Garnock-Jones 2257. О, P. Hebe albicans. —O. Shoot apex showing young lateral inflorescences and large apical 
vegetative bud, from Mt. Arthur, South Island, WELT 82555. —P. Shoot of plant from Cobb Valley. South Island, M. 
J. Bayly 5-31. —Q. Hebe pauciramosa, shoot of plant from Mt. Brewster, South Island, M. J. Bayly 1478. —W. Hebe 
speciosa, shoot from a plant cultivated at Otari-Wiltons Bush, Wellington, originally from Maunganui Bluff, North Island. 
—S. Hebe salicifolia, shoot of plant from Upper Wairau Valley, South Island, M. J. Bayly S-70. T-V. Hebe lycopodioides 
subsp. lycopodioides. —T. Habit of plant from same population as M. J. Bayly 1512, Mt. Nimrod, South Island. —U. 
Apical portion of vegetative branchlet, cultivated in Landcare Research Gardens, Lincoln, originally from Fish Lake, 
Tarndale, South Island, WELT 82551. —V. Terminal inflorescence of a plant from same population as M. J. Bayly 
771—773. Lake Tennyson, South Island. —W. Branchlets of Hebe armstrongii showing terminal infructescences, culti- 
vated in Landcare Research Gardens, Lincoln, originally from Nigger Stream, Canterbury, South Island, WELT 82552. 
—X. Small, light seeds of Hebe elliptica (typical of those of most Hebe), cultivated, Otari-Wiltons Bush, Wellington, 
North Island. Scale bars: A, B, = 1 em; С, G, I, L, N = 5 mm; E, V, W = 2 mm: J, K, U, X = 1 mm; О, P, Q, R, S 
= 2 cm. 








50 


of Chionohebe (e.g.. Fig. 5D) are united by several 
synapomorphies, including cushion habit with the de- 
cussate leaf pairs slightly offset to form а pseudo- 
spiral (Heads, 1994b), thick-walled bristle-like eglan- 
dular hairs, solitary flowers with erect long corolla 
tubes (e.g.. Fig. 5E). and corolla veins branching dis- 
tally in the tube (Garnock-Jones, 1993a). The Chio- 
nohebe A clade is sister to Parahebe planopetiolata, 
one of the species of Parahebe “Group B" of Ash- 
win (in Allan, 1961). With Parahebe planopetiolata 
they share а five-lobed corolla and hygrochastic 
capsule dehiscence, but the species in the Chio- 
nohebe B. clade also have these features. 

The Chionohebe B clade includes C. densifolia 
(Australian and New Zealand populations), plus 
Parahebe trifida from the informal Parahebe 
"Group B” of Ashwin (in Allan, 1961) (Figs. 3, 4). 
Parahebe birleyi is not included in this clade in all 
trees, and in the consensus tree forms a polytomy 
with it and several other clades in the complex. The 
species of the Chionohebe B clade are all alpine 
plants found in the southern region of the South 
Island. Parahebe trifida occurs in alpine flushes 
and snowbanks, / birleyi is a plant of nival rock 


ledges, and C. densifolia is found in a range of 


stony alpine habitats. Chionohebe densifolia is also 
found in the Kosciusco National Park in Australia. 
Chionohebe densifolia, P. trifida, and P. birleyi are 
all similar in appearance. Several of their shared 
characters are likely to be synapomorphies, includ- 
ing few-flowered inflorescences, large flowers with 
purple anthers, and presence of long glandular 
hairs on leaves. Other shared characters are also 
shared with the Chionohebe A clade, including old 
leaves withering and fading but retained on stems. 
All three species of the Chionohebe B clade are 
thought to form hybrids with cushion-forming spe- 
cies of Chionohebe A clade (Wagstaff & Garnock- 
Jones, 2000). There is little evidence from mor- 
phology to separate the Chionohebe A and B clades. 


Garnock-Jones (1993a), from a cladistic analysis of 


morphological and flavonoid data, proposed a more 
inclusive “Chionohebe” clade incorporating all the 


species of Chionohebe and Parahebe “Group B7 of 


Ashwin (in Allan, 1961). Such a grouping appears 
paraphyletic at best in this and earlier ITS studies 
(Wagstaff & Garnock-Jones, 2000). 

The Parahebe clade includes all the represen- 
tatives of the informal Parahebe “Group A" of Ash- 
win (in Allan, 1961) plus P. spathulata. Parahebe 
spathulata is anomalous in this clade, which oth- 
erwise has morphological support. from inflores- 
cence, floral, and flavonoid characters (Garnock- 
Jones, 1993a). A sister relationship between this 


clade and P canescens, the sole species of Ashwin's 
| 
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"Group C." has 9396 jackknife support. The infor- 
mal "Parahebe Groups А and C" share several flo- 
ral apomorphies (Garnock-Jones, 1993a). Parahebe 
brevistylis and P. linifolia (from Ashwin's “Group 
B^) form a weakly supported small clade, which is 
sister to the "Group А and C" species. The clade 
is represented in both the North and South Islands 
of New Zealand, where species occur їп well- 
drained soils associated with river banks, cliffs, and 
screes. Parahebe canescens is a creeping diminutive 
herb of South Island lake shores; its reduced fea- 
tures match convergent similarities seen in other 
plants associated with this habitat. The entire clade 
except for P. brevistylis and P. spathulata is united 
by floral features such as short corolla tubes, col- 
ored nectar guides (Fig. 5B), and stamen filaments 
narrowed at the base. In the case of P. brevistylis, 
the differences can be explained as losses of ad- 
aptations for insect pollination (Garnock-Jones, 
19765). Parahebe spathulata shares some features 
of habit and flower morphology with P. cheesemanii 
and might have an allopolyploid origin involving 
species from the Parahebe clade and the Р. chee- 
semanit lineage. 

The Heliohebe clade (Figs. 3. 4) was formerly 
recognized as Hebe "Paniculatae" in the informal 
classification of Moore (in Allan, 1901), and later 
segregated as a distinct genus by Garnock-Jones 
(1993b) (Table 1). It was also previously recognized 
as a distinct group in the key of Cheeseman (1925). 
Heliohebe includes five species that are found in 
northeastern parts of the South Island on rock out- 
crops, cliffs, and sometimes in grassland. Mono- 
phyly is well supported by two unique indels and 
several possible morphological apomorphies in- 
cluding an inflorescence that is a terminal, com- 
pound raceme or spike (Fig. 5C). protogyny (also 
evident in Fig. 5C, where styles are protruding from 
buds on the lowest inflorescence branches), sta- 
mens erect, anthers cream or yellow, seeds fusiform 
to irregular in shape and winged, and hemitropous 
ovules (Garnock-Jones, 1993b). 

The Hebe clade corresponds to Hebe sensu Moore 
(in Allan, 1961) with the exclusion of Heliohebe and 
members of the Leonohebe clade. The majority of 
species within the Hebe clade form a well-support- 
ed group (100% jackknife value), with weaker sup- 
port for Hebe macrantha (Fig. 5F, G) and H. petriei 
as sisters to this group (Figs. 3. 4). The Hebe clade 
is largely endemic to New Zealand, including many 
of its surrounding islands, with two species also 
extending to South America, and one species (Н. 
rapensis, not included in this analysis) endemic to 
Rapa Island (Fig. 6). The clade is large and mor- 
phologically diverse: unambiguous morphological 
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synapomorphies are difficult to identify, with ab- 
sences or reversals in some taxa. Wagstaff and Gar- 
nock-Jones (1998) suggest that synapomorphies 
may include: a shrubby or arborescent habit (е... 
Fig. 5M. R, S). large leaf bud (e.g.. Fig. 50), entire 
leaf margins, protandrous flowers, peltate placentas. 
acute capsule apices, 3/5 or 5/8 inflorescence phyl- 
lotaxis. Within the Hebe clade there is little reso- 
lution, but several relationships are worthy of note. 
Firstly. the "Connatae" of Moore [in Allan, 1901. 
represented by И. benthamii. H. petriei, H. epacri- 
dea, (Fig. 5H. |) and H. ramosissima| are ројурћу- 
letic, with some members closely related to “Bux- 
ifoliateae” (H. odora and H. pauciramosa). These 
relationships were first suggested by Wagstaff and 
Garnock-Jones (1998). and the addition of further 
taxa in this study continues to support this earlier 
assessment. Secondly, as found by Wagstaff and 
Wardle (1999), three of the cupressoid species, H. 
salicornioides, H. armstrongii (Fig. ЭМ), and H. an- 
nulata (which share possession of fused anterior 
calyx lobes and chromosome number of n = 21), 
form a well-supported clade. Thirdly, two of the 
Chatham Island endemics, H. chathamica and H. 
dieffenbachii, are sister species (Figs. 3, 4), with 


Map showing postulated dispersal of hebes from the main islands of New Zealand. 


the third Chatham Island endemic. H. barkeri, be- 
ing placed in a large polytomy that includes this 
grouping. Hebe macrantha, placed with weak jack- 
knife support at the base of the Hebe clade. lacks 
many of the previously mentioned synapomorphies 
for the group. This species has long held an isolated 
or ambiguous position, being placed by Moore (in 
Allan, 1961) in its own grouping, Hebe "Grandiflo- 
rae," and was included in Parahebe by Heads 
(1987, 1994b). For the present we suggest its re- 
tention in Hebe, which still leaves the genus, with 
the exclusion of Heliohebe (Garnock-Jones, 1993b) 
and the Leonohebe clade, monophyletic (see ‘Table 


1). 
TAXA OF UNCERTAIN AFFINITIES 


Many of the genera and subgeneric groupings 
historically recognized in the Hebe complex (e.g.. 
Heliohebe, cushion-forming species of Chionohebe, 
and some Hebe groups) are shown to be monophy- 
letic in this study. Significant exceptions are Par- 
ahebe “Group В” (of Ashwin in Allan. 1901) and 


species groups within Hebe [e.g.. sects. Hebe, Sub- 
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distichae, Glaucae, and the informal group “Con- 
natae" of Moore (in Allan, 1961)] (see Table 1). 
The species of Parahebe “Group B" are widely 
scattered, forming a grade among the New Zealand 
hebes, and few show any close relationship to other 
members within their informal group. Previous au- 
Allan, 1961; 


1993a) have suggested the Parahebe 


thors (Ashwin in Garnock-Jones, 
1976b, 
"Group B" species exhibit considerable morpholog- 
ical similarity, but now it seems this. similarity 
might derive from plesiomorphie character states 
(relative to the Hebe—Heliohebe clade), from con- 
vergent evolution of floral features associated with 
self-pollination (Garnock-Jones, 19705), and. per- 
haps from reticulate evolution in P spathulata (dis- 
cussed above). Garnock-Jones (1993a) suggested 
that Parahebe “Group В” and Chionohebe should 
be united, but this view receives no support from 
the ITS analyses (Wagstaff & Garnock-Jones, 2000, 
and herein). Wagstaff and Garnock-Jones (2000) in- 
ferred that the ancestors of Parahebe and Chionoh- 
ebe initially evolved in a montane or alpine envi- 
ronment, then subsequently radiated into lowland 
environments during episodes of Pleistocene gla- 
ciation, Extinction probably had a more profound 
effect on the basal lineages of Parahebe and Chio- 
nohebe, and this process further confounds our ef- 
forts to resolve relationships. 

The relationships among species of Hebe remain 
uncertain after ITS analyses (Figs. 3, 4). Branch 
lengths are too short for us to confidently derive an 
infrageneric classification from this study. | may 
be that speciation and diversification in the Hebe 
clade is too recent for ITS divergence to reveal its 


phylogenetic pattern. Reticulate evolution might 


also have clouded the molecular signal either 


through diploid hybrid speciation or allopolyploidy 
in this group where about 32% of the species are 
polyploid. 


ORIGIN, DIVERSIFICATION, AND DISPERSAL 


Two widely differing opinions have been pre- 
sented regarding the age and origins of the New 
Zealand hebes. Some authors (e.g., Skipworth, 
1973: Heads, 1994а) have proposed a Gondwanan 
origin to account for the present distribution of the 
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group, whereas others (e.g., Raven, 1973; Garnock- 
Jones, 1993a: Wagstaff & Garnock-Jones, 1998) 
have suggested the group has arrived more recently 
in the Southern Hemisphere and that dispersal has 
played an important role in shaping its distribution. 
The data presented here lend support to the latter 
proposition. 

A primary contribution of the present work is its 
assessment of relationships between Veronica and 
Australasian members of the Hebe complex. Al- 
though a close relationship between these two 
groups has long been assumed on morphological 
grounds (e.g.. Wettstein, 1891; Cheeseman, 1925; 
Raven, 1973), some authors (e.g., Hong, 1984) have 
directly opposed the notion that Australasian taxa 
are derived directly from within Veronica. Our anal- 
ysis of rbcL clearly supports a close relationship of 
Veronica to the Australasian genera, with the strict 
consensus including a well-supported (10096 jack- 
knife) Hebe—Veronica clade (Fig. 3). Analysis of ITS 
sequences shows the Australasian members to form 
a clade nested within a paraphyletic Veronica. This 
pattern of relationships is congruent with the notion 
that the Hebe complex is an Australasian radiation 
of Veronica. 

Our assessment of the time frame for the origin 
and radiation of the Australasian genera, and the 
New Zealand hebes in particular, relies on infer- 
ences from the fossil record, the geological and cli- 
matic history of New Zealand, and the distributions 
and ecological tolerances of extant species. For 
Hebe, the earliest appearance in the fossil record is 
in the Pliocene (Mildenhall, 1980), for Scrophular- 
iaceae it is in the mid Miocene (Tiffney, 1985), and 
for the whole of the Lamiales (sensu APG, 1998) it 
is in the mid Eocene (Muller, 1981). Although there 
is always the possibility that older fossils will be 
found, it would be inconsistent with this record to 
assume that divergence, either within the Austra- 
lasian Hebe complex, or between members of that 
group and Veronica, occurred in Gondwanan (Cre- 
taceous) or earlier times. 

We acknowledge that the divergence estimates 
presented in Figure 2 are crude and await further 
refinement. The paucity of the fossil record and the 
occurrence of undetected multiple substitutions on 





Figure 7. 


— 


Details of the natural elevational range of species of New Zealand hebes, overlaid on the strict consensus 


of trees produced by analysis of ITS sequences. Species are regarded as alpine if they occur predominantly in areas 


above the natural tree line. The other two elevational zones are less precisely defined with, on the two main islands of 


New Zealand, lowland corresponding to those areas below ea. 500 m above sea level; and montane being those between 
ca. 500 m and up to ca. 1000 m (but below natural tree line). Members of the basal grade of New Zealand hebes all 


occur in alpine or montane environments. 
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long branches are two potentially significant sourc- 
es of error, both of which could lead to inaccurate 
estimates of divergence times. Though crude, the 
estimates nonetheless provide intriguing compari- 
sons with patterns of ecological diversification and 
aspects of the geological history of New Zealand. 
The geological and climatic history of New Zea- 
land suggests that differentiation of the New Zea- 
land hebes is likely to have occurred in the late 
Tertiary. Members of the basal grade of New Zea- 
land hebes, in particular the Leonohebe clade, and 
the Chionohebe A and B clades (Fig. 3), are all 
alpine or montane plants, most occurring in areas 
above the natural tree line (Fig. 7). If the present 
ecological requirements of these groups are indic- 
ative of those of their past (i.e.. assuming that each 
lineage has not independently and recently adapted 
to alpine habitats, or that each has not seen selec- 
live extinction of lowland members), it can be in- 
ferred that early differentiation of New Zealand 
hebes occurred in alpine environments, with colo- 
nization of the lowlands being a secondary event. 
The evidence is that alpine environments have only 
existed in New Zealand since the Pliocene or latest 
Miocene, subsequent to the onset. of mountain 
building, in what was previously relatively low-ly- 
ing land (Flemming, 1979: Ollier, 1986). 
Prerequisite in any hypothesis of a late Tertiary 


origin for the New Zealand hebes is colonization of 


New Zealand by long-distance dispersal of ances- 
tral form(s). Assuming a minimum number of dis- 
persals or extinctions, the topology of cladograms 
derived from ITS sequences (Figs. 3. 4) suggests 
that. differentiation of New Zealand hebes followed 
a single colonization from either Australia or Eur- 
asia. This differentiation was succeeded by second- 
ary dispersal from New Zealand. 

As illustrated here and elsewhere (Wagstaff & 
Garnock-Jones, 1998, 2000) the morphology and 
ITS sequences (Fig. 3) of extant species with trans- 
oceanic distributions provide evidence of the ca- 
расцу of New Zealand hebes for long-distance dis- 
persal. Such transoceanic distributions are seen in 
four species, all of which are included in our ITS 
study group. and all of which are nested within the 


well-supported clade of New Zealand hebes. Of 


these species, Hebe elliptica and H. salicifolia nat- 
urally occur both in South America and southern 
New Zealand, whereas Chionohebe ciliolata and C. 
densifolia occur both in the South Island of New 


Zealand and southeastern Australia. Within each of 


these four species, populations separated by ocean 
gaps show no apparent morphological differentia- 
tion and only limited sequence divergence (Fig. 4). 
This suggests that the transoceanic disjunctions in 


the distributions of these species are relatively re- 
cent phenomena, and given that these species and 
all of their closest relatives occur within New Zea- 
land, that these disjunctions are products of long- 
distance dispersal of propagules from New Zealand. 

Godley (1967) suggested oceanic birds as likely 
vectors for the dispersal of seeds of H. elliptica (Fig. 
5X) and H. salicifolia from New Zealand to South 
America. Trans-Tasman dispersal of the two alpine 
Chionohebe species from New Zealand to Australia 
is less intuitively explained owing to the sexual di- 
morphism of one species (Delph, 1988, 1990), and 
their splash cup method of seed dispersal (Gar- 
nock-Jones, 1993a). The implied direction of dis- 
persal from New Zealand to Australia is also 
against the prevailing westerly winds but, as noted 
by Wardle (1978), weather conditions sometimes 
occur in which the usual direction of winds across 
the Tasman Sea is reversed. The occurrence of two 
independent. dispersals of Chionohebe from New 
Zealand to Australia may seem unlikely, but the 
alternative explanations are either an extended pe- 
riod of stasis in both morphology and ITS sequenc- 
es (assuming distributions produced by fragmenta- 
tion of Gondwana), or widespread extinction. in 
Australia (assuming dispersal in the opposite di- 
rection). 

Apart from the dispersal prerequisite to explain 
Iransoceanie species distributions, another six dis- 
persal events from the main islands of New Zealand 
are required to explain the current distribution of 
Hebe (Fig. 6). Most of the postulated dispersals are 
to New Zealand’s outlying islands, including one to 
the Pleistocene-age (Sykes, 1977) Кегтадес 15- 
lands (where Hebe breviracemosa is endemic), and 
three to islands of the New Zealand subantarctic 
(where H. elliptica and Н. odora have populations 
disjunct from those on the main islands of New 
Zealand, and H. benthamii is endemic). One dis- 
persal to the Chatham Islands has been postulated 
on morphological grounds (Moore in Allan, 1961; 
Garnock-Jones, 1976a: Wagstaff & Garnock-Jones, 
1998). This is partially supported here by analysis 
of ITS sequences (Fig. 4), which places two of the 
endemic species, H. chathamica and H. dieffenba- 
chii, as sister taxa, and the third, H. barkeri, in the 
polvtomy that includes the branch uniting the other 
two. A final dispersal, probably from the Chatham 
Islands (Garnock-Jones, 1976a, 1993a), is also pos- 
tulated to account for the distribution of H. rapensis 
(not included in our analysis), which is endemic on 
Rapa Island in French Polynesia. 

The presence of Parahebe in New Guinea is dif- 
ficult to explain. Here, as in the analysis of Wagstaff 
and Garnock-Jones (2000), ITS sequence data for 
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only one New Guinean species of Parahebe were 
included. That species, Р. vanderwateri, is nested 
within the New Zealand hebes (Fig. 3) with 99% 
jackknife support. The most parsimonious interpre- 
tation of the present data (assuming a minimum 
number of dispersals or extinctions) is long-dis- 
tance dispersal from New Zealand to New Guinea, 
as proposed by Wagstaff and Garnock-Jones (2000). 
Further sequence data for New Guinean Parahebe 
(of which 12 species are described) and Detzneria 
(monotypic and endemic) might provide a clearer 
picture of relationships between taxa from the two 


areas. 
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